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Abstract 

Background: The monocyte chemoattractant protein 1 (MCP-l) is involved in the recruitment of lymphocytes and 
monocytes and their migration to sites of injury and cellular immune reactions. In a Ghanaian tuberculosis (TB) 
case-control study group, associations of the MCPl -362C and the MCPl -2581G alleles with resistance to TB were 
recently described. The latter association was in contrast to genetic effects previously described in study groups 
originating from Mexico, Korea, Peru and Zambia. This inconsistency prompted us to further investigate the MCPl 
gene in order to determine causal variants or haplotypes genetically and functionally. 

Results: A 14 base-pair deletion in the first MCPl intron, intldel554-567, was strongly associated with protection 
against pulmonary TB (OR = 0.84, CI 0.77-0.92, Pcorrected = 0.00098). Compared to the wildtype combination, a 
haplotype comprising the -2581 G and -362C promoter variants and the intronic deletion conferred an even 
stronger protection than did the -362C variant alone (OR = 0.78, CI 0.69-0.87, Pnominai = 0.00002; adjusted Pgiobai = 
0.0028). In a luciferase reporter gene assay, a significant reduction of luciferase gene expression was observed in 
the two constructs carrying the MCPl mutations -2581 A or G plus the combination -362C and intldel554-567 
compared to the wildtype haplotype (P = 0.02 and P = 0.006). The associated variants, in particular the haplotypes 
composed of these latter variants, result in decreased MCP-l expression and a decreased risk of pulmonary TB. 

Conclusions: In addition to the results of the previous study of the Ghanaian TB case-control sample, we have 
now identified the haplotype combination -2581G/-362C/intldel554-567 that mediates considerably stronger 
protection than does the MCPl -362C allele alone (OR = 0.78, CI 0.69-0.87 vs OR = 0.83, CI 0.76-0.91). Our findings 
in both the genetic analysis and the reporter gene study further indicate a largely negligible role of the variant at 
position -2581 in the Ghanaian population studied. 



Background 

The monocyte chemoattractant protein 1 (MCP-l), also 
referred to as CCL2 (Chemokine [C-C motif] Ugand 2), 
is a member of the small inducible gene (SIG) family. 
CC-chemokines are characterized by two adjacent 
cysteine residues close to the amino terminus of the 
molecule. They are involved in the recruitment of lym- 
phocytes and monocytes and control migration of these 
cells to sites of cell injury and cellular immune reactions 
[1]. MCP-l is produced by different cell types in 
response to microbial stimuli [2]. MCP-l signals are 
transduced through the CCR2-receptor (chemokine [C- 
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C motif] receptor 2). Distinct microbial components are 
capable to induce expression of the CCR2 receptor and 
to initiate, dependent on the presence of MCP-l, target- 
oriented roaming of monocytes. 

The role of MCP-l in tuberculosis (TB) has been sub- 
ject of research since the early 1990 s. During the course 
of an infection with agents of the M. tuberculosis com- 
plex, MCP-l is predominantly produced by CD14+ 
blood monocytes and by distinct alveolar epithelial cells 
[3,4]. Elevated plasma MCP-l levels are found in TB 
patients [3], and the number of macrophages in bronch- 
oalveolar lavage fluids in eosinophilic pneumonia corre- 
lates with plasma MCP-l levels [5]. 

The gene encoding MCP-l {MCPl; MIM +158105) is 
located in the 17qll.2-ql2 chromosomal region. It con- 
sists of three exons and clusters with the loci CCL7, 
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CCLlh N0S2A, CCL3-S and CCR7. MCP-1 expression 
may be found in several conditions that are character- 
ized by infiltration of mononuclear cells. 

Genetic associations of MCPl gene variants with sus- 
ceptibility or protection against mycobacterial infection 
have been reported. Carriers of the MCPl -258 IG allele 
were at an increased risk of clinical TB in study groups 
from Mexico and Korea [4], Peru [6] and Zambia [7] 
compared to individuals carrying the alternative -2581 A 
allele. For MCPl -258 IG, several studies have demon- 
strated increased gene expression in vitro and elevated 
MCP-1 plasma levels in vivo [4,6,8-10]. In contrast to 
these observations, Japanese MCPl -258 lAA genotype 
carriers exhibited higher MCP-1 plasma levels than did 
carriers of the -2581 AG or GG genotypes [11]. While 
no effect of MCPl -2581 variants on TB susceptibility 
was observed in Brazilian, Chinese, Russian and South 
African study groups [12-15], MCP-1 -2581G was asso- 
ciated with protection against TB in a Ghanaian case- 
control study group [14]. Notably, the latter finding was 
in clear contradiction to the findings reported in Refs. 
[4,6,7], where an increased TB risk was attributed to 
MCPl -258 IG carriers. 

In the study of the Ghanaian TB case-control group, 
MCPl -258 IG was in weak linkage disequilibrium (LD) 
with another promoter variant, MCPl -362C (r^ = 0.27), 
which was even stronger associated with resistance to 
TB than MCPl -2581G. MCPl -362C has been shown 
to mediate increased transcriptional activity in vitro in a 
Caucasian study population [16]. Again, this finding is 
in contrast to the previous assumptions, namely that 
increased MCP-1 production might favour the occur- 
rence of TB [4,5,8]. 

These partly ambiguous findings prompted us to re- 
examine MCPl variants that might be involved in gene 
expression. According to the MCPl haplotype structure 
obtained from the innate immunity website https:// 
pharmgat.org/IIPGA2 eight variants that are located in 
the MCPl 5'-region, the first intron and in the 3'-UTR 
were selected and genotyped in our Ghanaian TB case- 
control group. Those genetic variants that showed the 
strongest associations with protection against TB were 
then subjected to a reporter luciferase gene assay in 
order to study gene expression. 

Results 

Alleles and genotypes 

Eight MCPl variants were genotyped in 2010 Ghanaian 
pulmonary TB cases and 2346 healthy control indivi- 
duals (Table 1). P values, including those of the previous 
study, were corrected according to the Bonferroni-Holm 
procedure [17] for the eight comparisons made. 

In Table 2, allelic associations of the eight newly geno- 
typed MCPl variants and of the MCPl -2581 and -362 



variants that were previously typed are given. The dele- 
tion of 14 bases length located in the first MCPl intron, 
intldel554-567, was associated with protection against 
pulmonary TB to a similar extent as were the promoter 
alleles -2581G and -362C (OR = 0.85, confidence inter- 
val [CI] 0.78-0.92, Pcorr = 0.00098, OR = 0.81, CI 0.73- 
0.91, Pcorr = 0.0012 and OR = 0.83, CI 0.76-0.90, P^orr = 
0.00015, respectively). 

The genotype frequencies did not deviate from 
Hardy- Weinberg equilibrium (HWE) among cases and 
controls. Trend tests were performed to compare the 
frequencies of genotypes of cases and controls in an 
additive model and results were adjusted for gender, 
age and ethnicity. The results are given in Table 3. As 
also observed in the computation of allelic associa- 
tions, intldel554-567 was in the trend test signifi- 
cantly associated with protection against TB (ORtrend 
= 0.84, CI 0.77-0.92, P^orr = 0.00098). In a genotype 
test where heterozygous and homozygous genotypes 
were individually compared to the wildtype, a strong 
association was seen for both heterozygous and homo- 
zygous carriers of the intldel554-567 deletion (OR = 
0.80, CI 0.70-0.91, Pcorr = 0.0063 and OR = 0.73, CI 
0.61-0.87, Pcorr = 0.0042, respectively), indicating a 
dominant genetic effect. The association of 
intldel554-567 was of similar strength as that in het- 
erozygous and homozygous MCPl -362C carriers in 
the previous study (Ref. [14]; ORtrend = 0.83, CI 0.76- 
0.91, Pcorr = 0.00017). Both variants were in strong LD 
(r^ = 0.82). intldel554-567 was also in weak LD (r^ = 
0.27) with the MCPl promoter variant at position 
-2581. Figure 1 shows the r^ values of pairwise LDs of 
all variants examined in the present and in the pre- 
vious study [14]. 

Stratification for mycobacterial species {M, tuberculo- 
sis vs. M. africanum) and phylogenetic lineages did not 
reveal any differences in the associations. Thus, possible 
confounding exerted by mycobacterial species or distinct 
genotypes was excluded. 

Haplotypes 

We focused on haplotypic combinations comprising the 
polymorphisms genotyped in our previous study, MCPl 
-2581A/G and -362C/G, and the deletion or wildtype 
(W) at intron 1 positions 554-567 (intldel554-567/W), 
because the variant alleles at these positions are in LD 
and associations of these variants are presented here 
and have been described previously [4,6,7,14]. As the 
combination -2581A/-362G/W occurred as the most 
frequent haplotype in our study population (frequency 
[/] = 0.55) it was referred to as wildtype reference in 
further comparisons (Table 4). 

The haplotype combinations -2581G/-362C/ 
intldel554-567 and -2581A/-362C/intldel554-567 were 
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Table 1 Variants selected for genotyping 



MCP1 


variant 


rs # 


Primer oligonucleotides 


Sensor/Anchor oligonucleotides 


-97569 


C/G 


rs9891330 


F-TCTGATGCAGACAGCGAG 
R-CACCTGGAGTCCCAG^ 


S-GCCTCCCCCACCCCCA 

A-ATAGCTGTCGGGAGAGTCTGTAmGAAAGAGAA 


-38600 




rs8075337 


f-otctgtgagcatoggt 
r-agaaacaaaaa™gggcatctac 


S-ACI ICI 1 1 IGCIGIGI 1 1 lAI 1 1 lAI 1 1 IC 
A-CCAACATCTGGAmCTOAGGGACAGmACA^ 


-1727 




rs39 17886 


F-TGGGATOTCCAGGAAACC 
R-ACCCAGCmCGTOGG 


S-GAAGAAGAGATACTGGAATGGAAACATCC 
A-GGGTGGGAGTCTCAGCACATCTACTATOTGTC 


lntl:96 


C/G 


rs41 507946 


F-TAAGGCCCCCTOTOTC 
R-CTGTGTGGTOGGCTCA 


S-CCCACAGTOTGCmAACGCTAC 
A-mCCAAGATAAGGTGACTCAGAAAAGGACAAGGG 


lntl:554-567 


W/D 


rs3917887 


F-TCCCCAGCTGATOTCC 
R-TGACTCAGmCCTATGCTGTA 


S-mAACCGCTCCTCOTC 

a-gtccgtc™atgacaotgtaggca™tctag 


+1542 


cn 


rs 13900 


F-GACCACCTGGACAAGCA 
R-AmCTOAGGCATAATGmCACATC 


S-IAGCI 1 ICCCCAGACACCCIGI IMA 
A-CACAACCCAAGAATCTGCAGCTAAC™^CCC 


+2413 


err 


rs39 17890 


F-ATGAGACCTGAAC™™mA 
R-TOACCCTAACATOAAATC 


S-GATCATOAGAAAGGAGAAGGAAGAGTGG 
A-AGCAAATACCTGGAGGTAGAAATGGTGATGATGTGTAC 


+2580 




rs41 343046 


F-GCCCACACCAATGTCAA 
R-CTGAATCTCTAAACATGGCAC 


S-AAGGGAmGAATGmGGGTGAAAAGATA 
A-ACTCAACTCTGTAGGTOAAAGGAAACGTOAGAA 



F, forward primer; R, reverse primer; S, sensor; A, anclior; Inti, intron 1; D, deletion; W, wildtype 



significantly associated with resistance to TB compared 
to the reference haplotype -2581A/-362G/W (OR = 
0.78, CI 0.69-0.87, P = 0.00002 and OR = 0.87, CI 0.78- 
0.96, P = 0.008, respectively; Table 4, Figure 1). The glo- 
bal P value, adjusted through 10 000 permutations, was 

Pglobal/adjusted = 0.0028. 



Reporter gene assay 

In order to test variant MCPl haplotypes with regard 
to their impact on gene expression, a luciferase repor- 
ter gene assay was performed. Figure 2 shows the plots 
of the Firefly Luciferase/Renilla Luciferase ratios (FL/ 
RL ratios) that were obtained for the constructs 



Table 2 Allelic associations 



MCPl 


allele 


cases n (frequency) 


controls n (frequency) 


OR 


CI 


p 

r nom 


p 

r corr 


-97569 


C 


1497 (0.38) 


1 745 (0.38) 


1 




0.99 






G 


2471 (0.62) 


2889 (0.62) 


1.00 


[0.92-1.09] 






-38600 


T 


1 730 (0.44) 


2107 (0.46) 


1 










C 


2228 (0.56) 


2463 (0.54) 


1.10 


[1.01-1.20] 


0.031 




-2581* 


A 


3256 (0.83) 


3692 (0.80) 


1 










G 


672 (0.17) 


932 (0.20) 


0.81 


[0.73-0.91] 


0.0002 


0.0012 


-1727 


A 


113 (0.03) 


149 (0.03) 


1 










T 


3525 (0.97) 


4133 (0.97) 


1.14 


[0.88-1.45] 


0.32 




-362* 


G 


2266 (0.58) 


2441 (0.53) 


1 










C 


1670 (0.42) 


2161 (0.47) 


0.83 


[0.76-0.90] 


0.000019 


0.00015 


lntl:96 


C 


967 (0.29) 


1092 (0.26) 


1 










G 


2401 (0.71) 


3128 (0.74) 


0.87 


[0.78-0.96] 


0.0055 




lntl:554-567 


W 


2432 (0.61) 


2646 (0.57) 


1 










D 


1586 (0.39) 


2028 (0.43) 


0.85 


[0.78-0.92] 


0.00014 


0.00098 


+ 1542 


T 


711 (0.18) 


933 (0.20) 


1 










C 


3291 (0.82) 


3745 (0.80) 


1.16 


[1.04-1.29] 


0.008 




+2413 


T 


134 (0.07) 


198 (0.07) 


1 










G 


1 750 (0.93) 


2652 (0.93) 


0.97 


[0.77-1.22] 


0.81 




+2580 


A 


397 (0.10) 


454 (0.10) 


1 










T 


3591 (0.90) 


4184 (0.90) 


0.97 


[0.84-1.12] 


0.71 





OR, odds ratio; CI, 95% confidence interval; P values are adjusted for age, gender and ethnicity; 

Pnom/ nominal P value; Pcom P value after Bonferroni-Holm correction; ^variants -2581 and -362 were originally genotyped and described in Ref. [15]; Inti, intron 



1 ; D, deletion; W, wildtype 
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Table 3 Genotype associations 



MCP1 


GT 


cases n (frequency) 


controls n (frequency) 


OR 


CI 


Pnom 


Pcorr 


ORtrend 


CI 


P P 

nom corr 


-97569 


CC 


278 (0.14) 


328 (0.14) 


1 








1.00 


[0.92-1.10] 


0.99 




CG 


941 (0.47) 


1089 (0.47) 


1.02 


[0.85-1.23] 


0.81 












GG 


765 (0.39) 


900 (0.39) 


1.01 


[0.84-1.22] 


0.92 










-38600 


CC 


629 (0.32) 


666 (0.29) 


1 








0.91 


[0.83-0.99] 


0.03 




CT 


970 (0.49) 


1131 (0.50) 


0.91 


[0.79-1.05] 


0.20 












TT 


380 (0.19) 


488 (0.21) 


0.83 


[0.70-0.99] 


0.03 










-2581* 


AA 


1355 (0.69) 


1472 (0.64) 


1 








0.81 


[0.73-0.91] 


0.0003 0.0018 




AG 


546 (0.28) 


748 (0.32) 


0.79 


[0.69-0.90] 


0.001 


0.006 










GG 


63 (0.03) 


92 (0.04) 


0.73 


[0.53-1.02] 


0.064 










-1727 


TT 


1 708 (0.94) 


1994 (0.93) 


1 








0.88 


[0.69-1.13] 


0.32 




AT 


1 09 (0.06) 


145 (0.07) 


0.87 


[0.67-1.12] 


0.27 












AA 


2 (<0.01) 


2 (<0.01) 


1.34 


[0.19-9.58] 


0.77 










-362* 


GG 


672 (0.34) 


654 (0.28) 


1 








0.83 


[0.76-0.91] 


0.000026 0.00021 




CG 


922 (0.47) 


1 1 33 (0.49) 


0.80 


[0.69-0.92] 


0.001 


0.008 










GG 


374 (0.19) 


514 (0.22) 


0.70 


[0.59-0.83] 


0.00005 


0.0004 








lntl:96 


CC 


145 (0.09) 


152 (0.07) 


1 








0.87 


[0.78-0.96] 


0.006 




CG 


677 (0.40) 


788 (0.37) 


0.90 


[0.70-1.16] 


0.4 












GG 


862 (0.51) 


1 170 (0.55) 


0.77 


[0.60-0.98] 


0.037 










lntl:554 


WW 


743 (0.37) 


734 (0.31) 


1 








0.84 


[0.77-0.92] 


0.00014 0.00098 


-567 


DW 


946 (0.47) 


1 1 78 (0.50) 


0.80 


[0.70-0.91] 


0.0009 


0.0063 










DD 


320 (0.16) 


425 (0.18) 


0.73 


[0.60-0.98] 


0.0006 


0.0042 








+1542 


CC 


1358 (0.68) 


1494 (0.64) 


1 








0.86 


[0.77-0.96] 


0.009 




CT 


575 (0.29) 


757 (0.32) 


0.84 


[0.73-0.95] 


0.007 












TT 


68 (0.03) 


88 (0.04) 


0.84 


[0.60-1.16] 


0.29 










+2413 


GG 


813 (0.86) 


1232 (0.86) 


1 








1.02 


[0.92-1.29] 


0.819 




GT 


124 (0.13) 


188 (0.13) 


0.99 


[0.78-1.28] 


0.99 












TT 


5 (<0.01) 


5 (<0.01) 


1.54 


[0.44-5-47] 


0.5 










+2580 


AA 


26 (0.01) 


21 (<0.01) 


1 








0.97 


[0.84-1.12] 


0.711 




AT 


345 (0.17) 


412 (0.18) 


0.69 


[0.38-1.25] 


0.22 












TT 


1623 (0.81) 


1886 (0.81 


0.70 


[0.39-1.25] 


0.23 











GT, genotype; Inti, intron 1; D, deletion; W, wildtype; OR, odds ratio; CI, 
nominal P value; Pcom P value after Bonferroni-Holm correction. ORtrend/ 
genotyped and described in [14]. 



95% confidence interval. P values are adjusted for age, gender and ethnicity. Pnom/ 
estimates of an additive genetic model; ^Variants -2581 and -362 were originally 



subjected to the assay. An overall ANOVA statistics 
revealed a significant difference between the FL/RL 
ratios (P = 0.0019). The calculated studentized range 
critical value in the post hoc pairwise comparisons for 
variable groups (Tukey-HSD test) was 4.03, and com- 
parisons of the construct carrying the -2581A/-362G/ 
W alleles with the -2581G/-362C/intldel554-567 and 
-2581A/-362C/intldel554-567 constructs yielded 
significant results that were above the studentized 
range critical value (Tukey-HSD test 4.53 and 5.44, 
respectively). 

The P values of a t-test that was calculated with the 
haplotype -2581A/-362G/W set as reference are given in 
Table 5. The constructs -2581G/-362C/intldel554-567 
and -2581A/-362C/intldel554-567 expressed the lucifer- 
ase gene to a significantly lower degree than did the 
wildtype construct -2581A/-362G/W (P = 0.02 and P = 
0.006, respectively). 



Discussion 

We have recently described an association of the MCPl 
-258 IG and -362C alleles with protection against TB in a 
Ghanaian study group [14]. The variants at these posi- 
tions were in weak linkage disequilibrium (LD; r^ = 0.27; 
Figure 1). In contrast to our observation of protection 
conferred by the MCPl -258 IG and -362C alleles, earlier 
research in study groups from Mexico, Korea, Peru and 
Zambia has attributed to MCPl -258 IG an increased risk 
to develop cUnical TB [4,6,7]. Moreover, MCPl -2581G 
and -362C were shown to enhance promoter activity in 
vitro in individuals from Korea {MCPl -2581G) and in 
Caucasians {MCPl -2581G and -362C), respectively 
[8,9,15], consistent with studies indicating that in pul- 
monary TB MCP-1 plasma levels are elevated [3,4]. 

Among the eight variants that were in LD with MCPl 
-362 and/or MCPl -2581 according to data available 
from NCBI and therefore subjected to genotyping in the 
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Figure 1 Pairwise linkage disequilibrium (LD) plots of MCP1 variants. Pairwise linkage disequilibrium (LD) plots of MCPl variants in the 
present and the previous study, rs numbers of variants that were newly genotyped are shaded. The haplotype consisting of variants -2581 A/G, 
-362C/G and intldel554-567/W are marked in turquoise blocks. 



present study we found the 14 base pair deletion in 
intron 1 (intldel554-567) associated with protection 
against TB similar to the ORs of MCPl -2581G and 
MCPl -362C. 

As the three variants are in weak LD {MCPl -2581 A/ 
G with MCPl -362C/G and MCPl intldel554-567/W) 
or in strong LD {MCPl -362C/G with MCPl 
intldel554-567/W) in the study group (Figure 1), haplo- 
types might explain more adequately than single muta- 
tions the genetic association and its relation to MCP-1 
production. All haplotypic combination comprising the 
variants MCPl -2581G/A, MCPl -362C/G and 
intldel554-567/W and occuring at frequencies (/) >0.01 
were considered (Table 4). The strongest association of 
protection against TB was with the haplotypic combina- 
tion -2581G/-362C/intldel554-567 (f = 0.19) compared 
to the wildtype haplotype -2581A/-362G/W {f = 0.54; 
OR = 0.78, CI 0.69-0.87, P = 0.00002). Carriers of the 
haplotypic combination -2581A/-362C/intldel554-567 



(f = 0.23) are slightly less, but still significantly protected 
against TB (OR = 0.87, CI 0.78-0.96, P = 0.008). 

To further examine the influence of haplotypes on the 
promoter activity, a reporter luciferase assay with con- 
structs comprising the -2581 A/G and -362C/G promo- 
ter variants, the first exon (149 bp) and the intronic 14 
bp deletion intldel554-567 or the wildtype sequence at 
these positions was performed. Only those combinations 
that occurred in frequencies >1% in the study popula- 
tion were included. The findings of the reporter gene 
assay corresponded to the results of the genetic analysis. 
A significant decrease of gene expression occurred in 
constructs carrying the -2581G/-362C/intldel554-567 
and -2581A/-362C/intldel554-567 combinations (P = 
0.02 and P = 0.006, respectively; Table 5). It may be 
inferred from the results of the reporter gene assay that 
both variants, MCPl -362C and intldel554-567, exert a 
reduction of the transcriptional activity, eventually 
resulting in lowered production of MCP-1. 
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Table 4 Associations of common haplotypes compared 
with a reference haplotype 



Table 5 Reporter gene assay 



Haplotype 




cases 


controls 


OR 


CI 


P 


-2582 -362 


Del 


n 

(frequency) 


n 

(frequency) 








A G 


W 


2281 
(0.57) 


2456 (0.52) 


1 






G C 


D 


677 
(0.17) 


928 (0.20) 


0.78 


[0.69- 
0.87] 


0.00002 


A C 


D 


881 
(0.22) 


1089 (0.23 


0.87 


[0.78- 
0.96] 


0.008 


A C 


W 


142 
(0.04) 


1 79 (0.04) 


0.84 


[0.67- 
1.07] 


0.15 



Del, intronic 14 bp deletion int1del554-567 OR, odds ratio; CI, 95% confidence 
interval. 

OR and P values refer to comparisons with the reference haplotype A/G/W. 
W, wildtype; D, deletion; global adjusted P value 0.0028 

Several mechanisms might be involved in the decrease 
of MCP-1 production. First, MCPl -362G constitutes a 
binding site for the Signal Transducer and Activator of 
Transcription 1 (Stat-1) [16]. Stat-1 is a transcription 
factor that enhances gene expression, and deprivation of 
Stat-1 binding through a loss of its binding site might 
reduce gene transcription, as seen in the luciferase 
assay. Another mechanism for a reduction of transcrip- 
tion is provided by the fact that Intronic deletions often 
cause a decrease of transcriptional activity [18,19] and 
influence RNA stability [20]. Lastly, the 14 bp deletion 
intldel554-567 in the first intron of the MCPl gene 
causes a loss of a predicted alternative splice site http:// 
zeus2.itb.cnr.it/ -'webgene/wwwspliceview_ex.html. While 
transcripts with several alternative splice sites appear to 
be more robust, a loss of transcription sites could pro- 
mote random degeneration in the nucleus [21]. It is. 



-2581A/-362G/WT -2581 G/-362C/DEL -2581 A/-362C/DEL -2581 A/-362C/WT 

Figure 2 Histogram illustrating the Firefly/Renilla (FL/RL) ratios. 

Histogram illustrating the Firefly/Renilla (FL/RL) ratios obtained after 
transfection of the four constructs; P < 0.05; P < 0.01; NS, not 
significant. 
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therefore, reasonable to ascribe a reduced MCP-1 pro- 
duction to the intldel554-567 deletion. 

Conclusions 

In addition to the results of the previous study of the 
Ghanaian TB case-control sample, we have now identi- 
fied the haplotype combination -2581G/-362C/ 
intldel554-567 that mediates stronger protection than 
does the MCPl -362C allele alone (OR = 0.78, CI 0.69- 
0.87 vs OR = 0.83, CI 0.76-0.91). Our findings in both 
the genetic analysis and the reporter gene study further 
indicate a largely negligible role of the variant at posi- 
tion -2581. 

The genetic risk of TB observed for variation in the 
MCPl promoter and in intron 1 is most likely conferred 
through an alteration of the MCPl expression, in line 
with the previous findings that increased MCP-1 pro- 
duction favours the occurrence of clinical TB. A similar 
observation was made in a mouse model for infections 
with both Listeria monocytogenes and M. tuberculosis, 
where increased MCPl expression in transgenic mice 
led to a 1 to 1.5 log greater sensitivity to infection [22]. 
It has been described that in MCPl deficient mice sub- 
jected to low-dose aerosol infection with M. tuberculosis 
Erdman the number of macrophages that enter the lung 
is decreased. As a consequence, these mice initially har- 
bour higher bacterial loads in their lungs compared to 
control animals, but eventually established a stable state 
of chronic disease [23]. No significant difference to 
MCPl wildtype mice in the susceptibility to intravenous 
infection was found [24]. It was also shown that mice 
that overexpress MCPl in their lungs exhibit increased 
uptake of M. tuberculosis BCG in dendritic cells com- 
pared to wildtype animals [25]. 

To date it remains unclear why high MCP-1 levels 
cause increased TB susceptibility in humans and how 
MCP-1 levels interact with the efficiency of the MCP-1 
gradient. Pertinent explanations could be that high sys- 
temic concentrations of MCP-1 would trigger the desen- 
sitization of receptors and reduce signal transduction or 
might lead to an adjustment and, as a consequence, to 
the neutralization of the chemoattractant gradient that 
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is required to escort sensitive monocytes to the sites of 
tissue damage. 

Methods 

Study group 

The study design and the enrollment procedure have 
been described in detail previously [14,26]. In brief, par- 
ticipants were recruited at the two major Ghanaian 
teaching hospitals in Accra and Kumasi and at addi- 
tional hospitals and polyclinics in these metropolitan 
areas and at regional district hospitals. 2010 HlV-nega- 
tive individuals with smear- and/or culture-positive pul- 
monary TB were recruited as cases. The control group 
consisted of 2346 individuals, from whom 1211 were 
unrelated personal contacts of cases and 1135 were 
community members from the adjacent neighbourhood 
or working contacts. The proportion of ethnic groups 
did not differ significantly between cases and controls. 
Participants belonged to the ethnic groups of Akan, Ga- 
Adangbe, Ewe and groups from northern Ghana, includ- 
ing Dagomba, Sissala, Gonja and Kusasi. The male-to- 
female ratio in the total study group was 1:0.58, and the 
mean age of participants was 33 years without gender 
differences. The characterization of phenotypes included 
documentation of the medical history of cases on struc- 
tured questionnaires, two independent examinations of 
non-induced sputum specimens, serological determina- 
tion of the HIV status, culturing and molecular differen- 
tiation of phylogenetic lineages of mycobacterial clades 
and posterior-anterior chest X-rays. Positive HIV test 
results were verified in an alternate test system. Fine- 
typing of genotypes by spoligotyping, /5'6110 fingerprint- 
ing and determination of drug resistance was performed 
as previously described [14]. TB-patients were included 
for specific treatment in the DOTS programme (Directly 
Observed Treatment Short-Course strategy) organized 
by the Ghanaian National Tuberculosis Programme. 

Of the control group, the medical history was obtained 
and a clinical examination was performed. Chest X-rays 
did not reveal any signs of actual or past pulmonary TB. 
In addition, a tuberculin skin test (TST, Tuberculin Test 
PPD Merieux, bioMerieux, Niirtingen, Germany) was 
performed. The TST was positive in 2217 controls and 
129 control individuals were TST-negative. 

Ethical approval of the study design was obtained by 
the Committee on Human Research, Publications and 
Ethics, College of Health Sciences, Kwame Nkrumah 
University of Science and Technology, Kumasi, Ghana, 
and the Ethics Committee of the Ghana Health Service, 
Accra, Ghana. Informed consent was given by study par- 
ticipants either by signature or, in case of illiteracy, by 
thumbprint in the presence of a witness. The aims of 
the study and the procedure of venous blood collection 
were explained before blood samples were taken. 



Variants selected for genotyping; genetic analyses 

According to the most recent data of the haplotype 
structure of MCPl obtained from the innate immunity 
database https://pharmgat.org/IIPGA2/PGAs/InnateIm- 
munity/CCL2/ we selected eight MCPl polymorphisms 
that are in LD with the MCPl -362 promoter variant 
which has previously shown the strongest association 
[14]. 

Table 1 lists the variants that were selected, including 
their rs numbers and PCR amplification primers as well 
as sensor/anchor nucleotides for LightTyper-based gen- 
otyping. Three variants are located in the promoter 
region, two in the first intron and three in the 3'-UTR. 

Standard methods were applied to extract DNA from 
full venous blood and genotypes of the MCPl variants 
were determined by fluorescence resonance energy 
transfer (LightTyper®; Roche Diagnostics, Mannheim, 
Germany) with dynamic allele specific hybridization. 

Databases and statistical analyses 

Demographic and self reported data was double entered 
into a Fourth Dimension database (San Jose, CA, USA). 
Genotype frequencies and odds ratios as well as Hardy- 
Weinberg equilibria (HWE) were calculated with the 
Stata 10 software (Stata Corporation, College Station, 
TX, USA) and logistic regression was applied to adjust 
for age, gender and ethnicities. Allelic and haplotype fre- 
quencies and associations were used to reconstruct hap- 
lotypes, calculated with the Unphased software (version 
3.1.4; http://www.mrc-bsu.cam.ac.uk/personal/frank/soft- 
ware/unphased). P values were adjusted through 10 000 
permutations. Haploview version 4.1 http://www.broad. 
mit.edu/mpg/haploview/ was used to calculate linkage 
disequilibria (LD, given as r^) and to generate the gra- 
phical output. The Tukey Honestly Significant Differ- 
ence test (Tukey-HSD test) was performed for post hoc 
comparisons of variable groups in the evaluation of the 
reporter gene assay. 

The power to determine a genetic effect (CaTS soft- 
ware; http://www.sph.umich.edu/csg/abecasis/CaTS/) 
with a genotype relative risk of 1.4 was, with 2010 cases 
and 2346 controls and assuming a disease allele fre- 
quency of 0.2, a prevalence 0.003 and a significance 
level of 1 X 10"^ was 89%. 

Reporter gene assay, engineering of constructs and 
transfection 

The PGL2-Control Vector (Promega, Mannheim, Ger- 
many) was used for cloning of all constructs of interest. 
Four fragments of the MCPl gene, each of 3569 bp 
length and containing the promoter, the first exon and 
the first intron, were PCR-amplified with primers 5'- 
caccaagaggagcttttcca-3' and 5'-gcgcacgcgtcctctgcactga- 
gatcttcct-3'. The MCPl -2581A/G and -362C/G variants 
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as well as the deletion intldel554-567/W were exam- 
ined. Only haplotype combinations occurring with fre- 
quencies >1% were subjected to the reporter gene assay. 
The following combinations were included: MCPl 
-2581A/-362G/W; -2581G/-362C/intldel554-567; 
-2581A/-362C/intldel554-567; -2581A/-362C/W. The 
Expand Long Template PGR System (Roche, Mannheim, 
Germany) was used for PCR-amplification. 

PGR conditions were: Initial denaturation (94°G, 2 
min), 10 amplification cycles (98°G, 10'; 60°G, 30'; 68°G, 
10"), 25 amplification cycles (98°G, 15'; 62°G, 30'; 68°G, 
20") and final elongation (78°G, 7"). SMAl and MLUl 
restriction sites at the 5' and 3'ends, respectively, were 
engineered on each PGR product. In an intermediate 
step, the fragments were gel-purified, ligated into a 
pGR-XL-TOPO plasmid (Invitrogen, Garlsbad, USA) 
and subsequently transfected into Top 10 cells (Invitro- 
gen, Garlsbad, USA) according to the manufacturer's 
instructions. After overnight incubation, cells were 
lysated and plasmids were digested with SMAl and 
MLUl restriction enzymes (New England Biolabs, Ips- 
wich, USA). The resulting fragments were then ligated 
into the PGL2 Gloning Vector and transfected into 
Top 10 cells. The final constructs were isolated using an 
EndoFree Plasmid Maxi Kit (Qiagen, Hilden, Germany). 

The Bio-RAD Gene Pulser Xcell system (Bio-Rad 
Laboratories Ltd., Hertfordshire, UK) was used for co- 
transfection of 6 x 10^ THPl-cells (German Resource 
Gentre for Biological Material, DSMZ [Deutsche Samm- 
lung fiir Mikroorganismen und Zellkulturen] , Braunsch- 
weig, Germany) with 0,5 \ig of the phRL-GMV vector 
and either 0,5 (ig of the pGL2-Gontrol vector or 0,5 \ig 
of one of the four plasmid constructs. Four hours after 
transfection, cells were harvested and luciferase activities 
were measured using a single tube Junior LB9509 
luminometer (Berthold Technologies, Bad Wildbad, 
Germany) and the Dual-Luciferase Reporter Assay Sys- 
tem (Promega, Mannheim, Germany). After a 10 second 
period of Firefly luminescence measurement, 100 ml Ix 
Stop & Glo Reagent that is supplied with the Dual-Luci- 
ferase Reporter Assay System kit were added and Renilla 
luminescence was detected in another 10 second mea- 
surement period. Ten independent transfections and 
measurements were performed for each construct. 
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